Abstract. Observations of infrared fine-structure lines 
provide direct information on the metallicity and ion- 
ization structure of Hn regions and indirectly on the 
hardness of the radiation field ionizing these nebulae. 
We have analyzed a sample of Galactic and Magellanic 
Cloud Hn regions observed by the Infrared Space Ob- 
servatory (ISO) to examine the interplay between stel- 
lar content, metallicity and the ionization structure of 
Hn regions. The observed [Siv] 10.5/[S in] 18.7 fim and 
[Ne in] 15.5/ [Neil] 12.8 /im line ratios are shown to be 
| highly correlated over more than two orders of magnitude. 

■ We have compared the observed line ratios to the results of 
j photoionization models using different stellar energy dis- 

■ tributions. The derived characteristics of the ionizing star 
£h \ depend critically on the adopted stellar model as well as 

, ^ ■ the (stellar) metallicity. We have compared the stellar ef- 
fective temperatures derived from these model studies for 
a few well-studied H II regions with published direct spec- 
troscopic determinations of the spectral type of the ion- 
^ ■ izing stars. This comparison supports our interpretation 
Tj" | that stellar and nebular metallicity influences the observed 

■ infrared ionic line ratios. We can explain the observed in- 
| crease in degree of ionization, as traced by the [S iv]/[S in] 
. and [Ne in] /[Neil] line ratios, by the hardening of the ra- 
j diation field due to the decrease of metallicity. The impli- 

f^*) • cations of our results for the determination of the ages of 
| starbursts in starburst galaxies are assessed. 
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1. Introduction 



Milky Way by Martin-Hernandez et al. (2002[) and in the 
Large and Small Magellanic Clouds by Vermeij & van dei 



and nebular ineLalhciL^ are iiiLunaLely inLei woven in con 

trolling the ionization structure of Hii regions. Hotter 
stars will ionize trace elements to higher ionization stages. 
Metallicity, on the other hand, modifies the spectral en- 
ergy distribution (SED) of the ionizing stars and thus, 
influences the ionization structure indirectly. In partic- 
ular, line blocking/blanketing is directly coupled to the 
metallicity and controls the stellar wind of massive stars. 
The stellar wind, in its turn, also modifies the SED of the 
star. The interlocking of these different effects hampers 
the determination of the stellar content from ionic line 



The s^LLicd Lvpe of fcfae ionizing s Lais and Lhe sLellai ggg (2QQ ^ haye nQw provided ug wit h a 'natural' sam 



pie of H II regions over a wide range in metallicity. Such 
a sample allows us to study the effect of stellar type and 
metallicity on the ionization structure. Moreover, while 
the ionizing stars in most H n regions are heavily extincted 
in the visible by tens to hundreds of magnitudes, the open- 
ing up of the infrared window has favoured in recent years 
the use of infrared spectra to infer the spectral types of 
the ionizing stars. An important effort has now started to 
identify the ionizing stars of Galactic Hii regions in the 
K-band (Hanson et al., 1996; Watson & Hanson, 1997) 



observations. This has clear repercussions oil Lhe study 
of starbursL regions hi galaxies as well, where Lhe sLellar 

content is used as an indicator of the starburst age. Cur- 
rent models of stellar atmospheres have s tarted to take 
the effects o f metallicity into accoun t (e.g. Schaerer fc de 
Koter, 1997| ; Pauldrach et al., 2001 ). Such models can fit 



and the first results are now becoming available (Kapei 
et al., 2002a| ,[b|) 



the observed line rat i os in H H regions reaso nable well (e.g 



Giveon et al., 200 1| ; [Morisset et al., 2002j ). The derived 



spectral types for the ionizing sources, however, depend 
on the stellar model adopted and can be very different. 

The Infrared Space Observatory (ISO, |Kessler et al 
1996| ) has provided a powerful tool for the study 
of the ionization structure of H n regions (jGivcon 



|et ah, 2001 ; Martin-Hernandez et al., 2002 , hereafter 



MHP02). The ratios of fine-structure lines originating in 
the ionized gas (for instance, [Arm] 9.0/[Arii] 7.0 /xm, 
[Siv] 10.5/[S in] 18.7 /xm and [Nem] 15.5/[Neii] 12.8 /xm) 
are widely used to constrain the properties of the ionizing 
stars in these nebula (e.g. Watarai et al., 199^ ; Takahashi 
et al., 2000; Morisset et al., 2002| ). The systematic studies 
of the ionization structure of H n regions via these ionic 
line ratios as a function of Galactic center distance in the 
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In this paper, we examine this interplay between stellar 
content, metallicity and ionization structure of H II regions 
based upon the ISO sample. This paper is structured as 
follows. Sect. H describes the combined sample of Galactic 
and LMC/SMC Hii regions; Sect. ^ shows the correlation 
between the observed [Nem]/[Neii] and [Siv]/[Siii] line 
ratios. Sect. || and Sect. || investigate the influence of the 
SEDs and the stellar/nebular metallicity, respectively, on 
the ionization structure. Sect. ^ discusses the implications 
for the interpretation of the spectra of starburst galaxies. 
Finally, Sect. [?] discusses and summarizes the results. 

2. The sample 

We combined the ISO Short Wavelength Spectrometer (|d| 
jGraauw et al., 1996| ) observations of the Galactic and Mag- 
ellanic Cloud H II regions as described in the catalogues by 
Peeters et al. (2002| ) and |Vcrmcij et al. (2002| ). These cat- 
alogues present the atomic fine-structure line fluxes of 43 
Galactic nebulae and 12 Magellanic Cloud sources. The 
first analysis of the observed fine-structure lines is pre- 
sented in MHP02 and |Vermeij fc van der Hulst (2002|) . 



This analysis includes, among others, the derivation of el- 
emental abundances and a discussion on the ionization 
state of the nebulae. 

From this combined set of Hn regions, we selected 
those for which both information on their elemental abun- 
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Table 1. Selected Hn regions in the Milky Way, Small 
Magellanic Cloud (SMC) and Large Magellanic Cloud 
(LMC). Given are the object name, the distance in kpc 
of the Galactic objects from the center of the Galaxy and 
the coordinates of the ISO pointing. 



Object 



Rgal RA (J2000.0) S (J2000.0) 



(kpc) 



(h ,m ,s) 



r, ") 



Milky Way: 

IR 01045+6506 
IR 02219+6125 
IR 10589-6034 
IR 11143-6113 
IR 12063-6259 
IR 12073-6233 
IR 12331—6134 



IR 15384—5348 
IR 15502-5302 
IR 17221-3619 
IR 17455-2800 
IR 18116-1646 
IR 18317-0757 
IR 18434-0242 
IR 18479-0005 
IR 18502+0051 
IR 19598+3324 
Dr21 

IR 21190+5140 
IR 23030+5958 



13.8 
11.0 
9.5 
9.7 
9.3 
10.1 

~inr~ 

4.6 

5.2 

0.5 

4.3 

4.5 

4.6 

7.5 

4.7 

9.8 

8.6 
12.7 
11.4 



01 07 50.7 

02 25 44.6 
11 00 59.8 

11 16 33.8 

12 09 01.1 
12 10 00.3 
VI 36 01.9 
lb 42 IV. 1 
15 54 06.0 
17 25 31.7 

17 48 41.5 

18 14 35.2 
18 34 24.9 
18 46 04.0 
18 50 30.8 
18 52 50.2 
20 01 45.6 

20 39 00.9 

21 20 44.9 
23 05 10.6 



+65 21 21.7 
+62 06 11.3 
-60 50 27.1 
-61 29 59.4 
-63 15 54.7 
-62 49 56.5 
-61 51 03.9 
-53 58 31.5 
-53 11 36.4 
-36 21 53.5 
-28 01 38.3 
-16 45 20.6 
-07 54 47.9 
-02 39 20.5 
-00 01 59.4 
+00 55 27.6 
+33 32 43.7 
+42 19 41.9 
+51 53 26.5 
+60 14 40.6 



SMC: 

N66 
N81 



00 59 03.7 

01 09 13.6 



-72 10 39.9 
-73 11 41.1 



LMC: 






N4A 


04 52 08.4 


-66 55 23.4 


N83B 


04 54 25.2 


-69 10 59.8 


N11A 


04 57 16.2 


-66 23 18.3 


30Dor#l 


05 38 33.5 


-69 06 27.1 


30 Dor#2 


05 38 35.5 


-69 05 41.2 


30Dor#3 


05 38 46.0 


-69 05 07.9 


30 Dor#4 


05 38 54.2 


-69 05 15.3 


N160A1 


05 39 43.3 


-69 38 51.4 


N160A2 


05 39 46.1 


-69 38 36.6 


N159-5 


05 40 02.4 


-69 44 33.4 



higher precision than other elements because it is less 
sensitive to ionization correction factors and assumptions 
regarding electron densities (see MHP02 for a complete 
discussion). Moreover, as a primary element, the abun- 
dance of neon closely follows that of oxygen (e.g. Henry 
fc Worthey, 1999| , and references therein). As far as the 



choice of the ionization tracers is concerned, the choice of 
the [Siv]/[Siii] and [Nein]/[Neii] line ratios was an obvi- 
ous one as [Arm]/ [Am] has been detected in only a few 
LMC Hn regions. 

In total, our sample consists of 20 Galactic Hn re- 
gions located at Galactocentric distances between and 
15 kpc, 7 Hn regions in the LMC (including 4 posi- 
tions in 30Doradus) and 2 in the SMC. Also included are 
the SWS observations of SgrA*, our own Galactic Cen- 
ter ( Lutz et al., 1996| ), and the Orion nebula, for which 



an "integrated" infrared spectrum was obtained by 3imp 



dances (through Ne/H) and ionization state (through the 
[S IV] 10.5/ [S III] 18.7 /xm and [Ne III] 15.5/[Nell] 12.8 pm 
line ratios) could be derived. The selected Hn regions 
are listed in Table [l], where the source name, the dis- 
tance of the Galactic objects from the center of the Galaxy 
and the ISO pointing coordinates are given. We chose to 
use Ne/H as a tracer of the nebular metallicity. Ne/H 
can be determined from the IR fine-structure lines to 



son et al. (1998 ) using the Midcoursc Space Experiment 
(MSX) satellite. The MSX field of view (6'x 9') included 
most of the Orion nebula. 



3. Variations in the ionization structure 

The ionization structure of a nebula depends basically on 
the shape of the SED and the nebular geometry. In the 
case of an ionized sphere of constant gas density and fill- 
ing factor, the geometrical effect can be defined by the 
ionization parameter U = Qh/ (4irR 2 nc), where Qh is the 
number of stellar photons above 13.6 eV emitted per sec- 
ond, R is the Stromgren radius, n is the gas density and 
c is the speed of light. The ionization structure can be 
traced by the ratio of successive stages of ionization X +l 
and X +l+1 of a given element. Such a ratio depends, for 
a given U , on the number of photons able to ionize X +l 



flVilchcz fc Pagel, 1988| ). 

in Fig. |lj, we present the relation between the 
[Siv]/[Siii] and the [Nein]/[Neii] line ratios for the com- 
bined sample of Hn regions in the Milky Way and the 
LMC/SMC. The line ratios have not been corrected for 
extinction. Because the lines [Neil] 12.8 and [Nem] 15.5 
/j,m are very close in wavelength, the [Ne in] / [Ne n] is prac- 
tically insensitive to extinction. However, the differential 
extinction between the [Siv] 10.5 and [S in] 18.7 /j,m lines 
(~ 0.08Ak if a standard extinction curve is considered, 
see MHP02) can affect the [Siv]/[Siii] line ratio. The ef- 
fect of this differential extinction is indicated in the lower 
right corner of Fig. [I] for a typical Ak of 2 magnitudes 
(MHP02). 

The data shown in Fig. [I] span a range in ionization of 
more than 2 orders of magnitude. The LMC/SMC points 
nicely overlap the Galactic trend at the high ionization 
end. The least squares fit to the data, displayed as a dot- 
ted line, gives a slope of 1.04 ± 0.05 and an intercept of 
0.61 + 0.04. Fig. [l| also illustrates the change in the degree 
of ionization of the Galactic H II regions with increasing 
distance from the Galactic center. The more highly ion- 
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Fig. 1. Relationship between the 

[Siv] 10.5/ [S in] 18.7 /^m and [Nelii] 15.5/[Neil] 12.8 /im 
line ratios for the combined sample of Hn regions. 
Indicated by various symbols are the Galactic nebulae 
at i?Gai< 7 kpc (solid triangles), the Galactic nebulae 
at i?Gai> 7 kpc (solid circles), the LMC nebulae (open 
triangles, except 30Doradus, which is indicated by an 
open square) and the SMC nebulae (open diamonds). To 
avoid confusion, only the averaged location of the 4 SWS 
positions in 30Doradus is indicated. The positions of 
Sgr A* ( Lutz ct al., 1996| ) and the Orion nebula ( [Simpson 



et al., 1998) are indicated by a solid star and a solid 



square, respectively. The dotted line is a least squares 
fit to the data. A typical error bar of the ISO data is 
given in the upper left corner. The arrow in the lower 
right corner indicates the extinction correction due to a 
Ak = 2 mag. 



ized objects are located at larger Galactocentric distances, 
while the less highly ionized objects, including Sgr A*, are 
located at smaller distances. In line with the known metal- 



licity gradient in the Galaxy (e.g. see review by Henry 
& Woithey, 1999; Rolleston et al., 2000, and references 



therein), the location of the high metallicity Galactic Cen- 
ter at the low ionization end of the correlation and the low 
metallicity LMC/SMC Hn regions at the high ionization 
end suggests that metallicity is somehow involved in the 
changes of ionization degree, and hence, in the hardening 
of the radiation field, shown in Fig. |l|. 

We consider in the next sections both the influence of 
the SED and the metallicity on the ionization structure of 
H II regions. 



4. Influence of the SEDs 

In order to investigate the dependence of the [Ne iil]/[Ne n] 
and [Siv]/[Siii] line ratios on the SED, a set of nebu- 
lar models has been calculated with the photoionization 
code CLOUDY ( |Fcrland ct al., 1998j ) using MICE, the 
IDL interface for CLOUDY created by H. Spoorf]. The 
parameters varied in the modelling are the stellar effective 
temperature {T c g) and the ionization parameter ([/). We 
compute nebular models for a static, spherically symmet- 
ric, homogeneous gas distribution with one ionizing star in 



the center. An inner cavity with a radius equal to 10 cm 
is set, while the outer radius of the H n region is defined by 
the position where the electron temperature, T c , reaches 
100 K. Two different stellar atmosphere models are taken 
to describe the SED: the CoStar models by [Schaerer fc de 
Koto CI 997ft and the models hv |Pa,nldra,ch ct al. fjgOjE 
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Fig. 2. Dependence of the [Siv]/[Siii] and [Nem]/[Neii] 
line ratios on the stellar T e s (crosses) as given by the 
CoStar (solid lines) and Pauldrach (dashed lines) SEDs. 
Only photoionization models for nn=10 3 cm -3 are plot- 
ted. The observed line ratios for reference sources (Sgr A*, 
G29.96-0.02, the Orion nebula, IRAS 12073-6233, N66 
and 30 Dor) are indicated by solid circles. The dotted line 
is the least squares fit to the data plotted in Fig |l|. 



Previous studies of these ionization ratios have concen- 



trate d on either the CoStar models (e.g. Morisset et al 



http : / / www. astro . rug. nl /~spoon/ mice .html 



2001 



2002 1, MHP02)) or on the Pauldrach models ( |Givcon ct al 



). For a detailed study of the influence of these and 
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other stellar SEDs on the ionization ratios, we refer to the 
work by Morisset et al. (in prep.). We use stellar models 
for main sequence (dwarf) stars, i.e. CoStar models A2, 
B2, C2, D2, E2 and F2, and Pauldrach models D-35, D-40, 
D-45, D-50 and D-55. The metallicity of the star and the 
nebula is chosen to be solar. The number of hydrogen ion- 
izing photons emitted by the central source is fixed to the 
typical value of Qh = 10 48 photons s _1 . The parameter U 
is varied by changing the hydrogen density («h). For each 
of the stellar atmosphere models selected, which span a 
range in T g from approximately 30 to 55 kK, photoioniza- 
tion models are calculated for n H = 10 2 , 5 x 10 2 , 10 3 , 5 x 10 3 
and 10 4 cm- 3 . 



Several points can be gleaned from inspection of Fig. |[ 
(a) The [Nem]/[Neii] and [Siv]/[Siii] line ratios 
strongly depend on the stellar T e g, but they can no longer 
discriminate between models when most of the Ne and S is 



in the higher ionization stage (see also Givcon et al., 2001 




o 



20 



CoStar: 36000 K 
CoStar: 46000 K 
Pauldrach: 45000 K 



MHP02). When this happens, these line ratios remain con- 
stant even with large changes in the SED. When using the 
CoStar SEDs, both [Nem]/[Neii] and [Siv]/[Siii] become 
insensitive to the stellar T e g for T e g> 42 kK, while when 
using the Pauldrach SEDs, the hardness of the stellar ion- 
izing flux changes less rapidly with T e g, and it is not until 
T e g> 50 kK that the star fully ionizes S ++ and Ne + . 

(b) The slope of the observed log([Neiil]/[Neii]) ver- 
sus log([Siv]/[Sm]) is practically equal to 1 (cf., Fig |l]), 
suggesting that the spectral hardening affects equally the 
range of ionizing energies for S in (hv > 35 eV) and Ne n 
(hv > 41 eV). The model trends are, however, slightly 
steeper than the data, i.e. for a given T e g the models 
overpredict [Neni]/[Nen] (underpredict [S iv]/[S in]) by 
a factor up to 5 for a giv en |"S iv]/|"S in] ([ Ne Hl] /|Ne ill). 
As has been discussed by |Ocy et al. (2000| ) and Morisset 
ct al. (2002[) , there is an indication that the CoStar models 
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Fig. 3. The number of ionizing photons in cm "s 
a certain energy from the CoStar B2 (36000 K), D2 (46000 
K) and Pauldrach D-45 (45000 K) SED models is plotted 
as a function of this energy. The curves are normalized so 
that the total number of hydrogen ionizing photons, Qh, 
is 10 48 photons s — 1 , the value used in the photoionization 
models. These curves define, for a given ionization param- 
eter U, the ionization structure in the nebula. Labeled on 
top of the figure are the ionization energies for H°, Ne°, 
S+, S++, Ne+ and He+. 



The output of this grid of models for nn=10 3 cm -3 
is presented in Fig. |[ The variation of [Ne ill] /[Neil] and 
[Siv]/[Siii] with T c ff is shown separately for the two se- 
lected SED models. An increase in density of two orders 
of magnitude produces an increase of the line ratios of 
at most 0.5 dex approximately along the lines drawn in 
Fig. ||. For reference, the observed line ratios of the Galac- 
tic sources SgrA*, G29.96~0.02 (IRAS 18434-0242), the 
Orion nebula and IRAS 12073-6233, the SMC Hn region 
N66, and 30Doradus are indicated by solid circles. SgrA* 
and IRAS 12073 delimit the observed parameter space for 
[Ne III] /[Neil] and [Siv]/[Siii] as they are, respectively, 
the lowest and highest ionized source in Fig. |l|. 



overpredict the stellar spectrum at these ionizing energies. 

(c) When using the CoStar SEDs to describe the emer- 
gent flux, practically the entire observed ionization range 
spanning more than 2 orders of magnitude can be mod- 
eled by a rather narrow T e g range between ~ 33 kK and 
~ 42 kK. This effect is less dramatic with the Pauldrach 
SEDs, for which the entire range up to 55 kK is suitable 
for modelling the observed line ratios. 

We see that large differences appear when we compare 
the predictions by the CoStar models to those by the Paul- 
drach ones for a given T c ff. These differences are primarily 
caused by the different treatment of the line blocking and 
blanketing between the two SED models. This causes the 
Pauldrach models to have a much softer spectrum than the 
analogous CoStar star at a similar T c g. Note, for instance, 
the case of the CoStar SED with T Gfr =36000 K (model B2) 
and the Pauldrach SED with T off =45000 K (model D-45) 
in Fig. |^, which predict practically the same values for 
[Siv]/[Sm] and [Ne in] /[Neil]. 

It is instructive to calculate the integrated ionizing 
photon flux above a certain energy as a function of this 
energy for these two models (CoStar B2 and Pauldrach 
D-45). The number of ionizing photons in cm _2 s _1 with 
energies larger than hvi is defined as q Vi ~ J u F v / (hv)dv, 
with F v being the astrophysical flux in erg s _1 cm~ 2 Hz -1 . 
The results are shown in Fig. |^. For a significant compari- 
son with the photoionization model results, the total num- 
ber of hydrogen ionizing photons has been scaled to 10 48 
photons s _1 , the value used in the models. We note that 
models differ in their stellar radii. The [Siv]/[Siii] and 
[Ne m]/[Ne n] line ratios are basically sensitive to the slope 
of the curve between the energies able to produce and ion- 
ize S++ and Ne+, i.e. between 23.3 and 34.8 eV, and 21.6 
and 41.0 eV, respectively. The slopes of the CoStar B2 and 
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Table 2. Comparison of the stellar properties predicted via the [Nem]/[Neii] and [Siv]/[Siii] line ratios with the 
spectral type derived directly from infrared/optical stellar spectra. Given are the object name, the observed spectral 
type of the ionizing star(s) with its estimated T e s, the predicted stellar properties by using either the CoStar or 
the Pauldrach SEDs and the reference of the observed spectral type. The objects are sorted by increasing degree of 
ionization. 



Observed T off t 



Predicted T cff t 



Object name 


Spectral Type 


(a) 


(b) 


CoStar 


Pauldrach 


Ref. 


G 29.96-0.02 


05-03 


46-51 


42-49 


33-36* 


35-40 


1 


Orion nebula 


06 


44 


39 


33-36 


40-45 


2 


N66 


06.5 and earlier 


> 42 


> 37 


33-42 


>40 


3 


30 Doradus 


04-03 


49-51 


43-49 


36-42 


>45 


4,5 


IRAS 12073-6233 


07.5-06 


40-44 


36-39 


> 36 


>45 


1 



(a) T e ff derived from the observed spe ctral type using the calibration by Vacca et al. (1996); (6) T e g derived from the observed 



spectra l typ e using the calibration by Martins et al. (2002); (f) T e s in kK; (*) The detailed photoionization model by 



et al. (fo02j), which uses the CoStar SEDs, p redicts a T cff of ~ 30 l \ kK. REFERENCES: 



communication; (2) SIMBAD database; (3) Massey et al. (1989); (4) Walborn fe Blades (1997); (5) Massey fc Hunter (199E 



Morisset 



Kapcr et al. ( 2002a ^) and priva te 



the Pauldrach D-45 SEDs are very similar between these 
energi es; hence, their match in the predicted [S iv]/[S in] 



and [N fciiil/[Nc ill line ratios. The ionization structure of a 
nebula with such values of [Siv]/[Siii] and [Ne m]/[Ne n] 
can, therefore, be equally well described by any of these 
two SEDs. However, the interpretation of this ionization 
structure in terms of the stellar properties, e.g. T e ff, de- 
pends on the stellar model adopted. For comparison, we 
also show the shape of the hotter CoStar model D2, with 
a T cff of 46000 K. The ionizing flux of this hotter CoStar 
model starts to be higher than the corresponding Paul- 
drach model above ~ 30 eV. This hotter CoStar model 
will produce, consequently, a larger amount of S +3 and 
Nc++. 

In this respect, closely related to the problem of de- 
termining the T ff of the ionizing sources is that of their 
spectral type classification. It is illustrative to consider 
at this point the cases of G29.96— 0.02, the Orion neb- 
ula, IRAS 12073-6233, N66 and 30 Doradus, for which the 
spectral types of the ionizing star are known. The Orion 
nebula is ionized by the Trapezium cluster, the dominant 
source of ionizing radiation being the 06 star 1 COri. 
30 Dor is known to be powered by a stellar cluster con- 



ing star in G29 is confirmed to be 03-05 (we n ote that 
lower resolution K-band spectral observations by Watson 



taming a large number of 03-04 stars (Walborn & Blades 



1997; Massey & Hunter, 1998). N66 is the most prominent 
Hn region in the SMC. It contains a rich stellar cluster of 
at least 33 O stars of which 11 are of type 06.5 or earlier 



(Massey et al., 1989). Recent high resolution K-band spec- 



tral observations of G29 and IRAS 12073 ( [Kaper et al 
2002a|]p[ private communication) allow the determination 
of the spectral type of the ionizing stars based mainly 
upon the presence and strength of photospherical lines of 
ClV and Nlli (Hanson et al., 1996). These observations 



indicat e that IRAS 12073 is ionized by a cluster of young 
06-Q7 I.5 stars, while the spectral type of the main ioniz- 



fc Hanson 1997 classified this star as being 05-08). 

As mentioned above, the interpretation of the 
[Siv]/[Siii] and [Ne m]/[Ne n] line ratios observed for 
these objects in terms of the stellar T c g depends on the 
stellar model used to describe the ionization structure 
traced by these ratios. The subsequent spectral classifi- 
cation of the ionizing source depends on the spectral type 
attributed to the SED used in the modelling of the H n 
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Fig. 4. The [Nein]/[Nen] line ratio is compared to the 
observed spectral type for G29.96— 0.02, the Orion nebula, 
IRAS 12073-6233, N66 and 30 Doradus. The numbers in 
brackets indicate the Ne/H abundance by number in units 
of 10 -4 . References for the Ne/H abund ances are: G29 and 
IRAS 12073 (MHP02 ); N66 and 30 Dor flVermeij fc van dei 
Hulst, 2002| ): Orion ( |5impson et al., 1998| ). 
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Fig. 5. Relationship between the [Ne in] /[Neil] line ra- 
tio and the metallicity indicator Ne/H for the combined 
sample of Hn regions. Symbols are as in Fig. [fl. The data 
collected on starburst galaxies (open stars) are compared 
to the Galactic and LMC/SMC Hn regions. The paral- 
lel lines delineate the parameter space span up by the 
data. A typical error bar for log([Ne ill] /[Neil]) is shown 
in the left bottom corner. No error bar for Ne/H is indi- 
cated as the main uncertainty (which can be as large as 
a factor of 2 for the Galactic sources, see Sect. |) comes 
from the abundance derivation method itself and is very 
hard to determine. The solid lines show the results of the 



photoionization models calculated by Givcon ct al. (2001 ) 
using Pauldrach SEDs for subsolar, solar and supersolar 
stellar and nebular metallicities. 



region. We can try to compare the T e g determined from 
the observed line ratios by using either CoStar or Paul- 
drach SEDS (cf., Fig ||) with the observed spectral type 
of the ionizing star(s). To do so, we need a T e ff-spectral 
type calibration for O stars. The T e s calibration of the dif- 
ferent spectral types is based on modelling the observed 
photosphcric hydrogen, helium and metal lines and it is by 
no means certain. For instance, the commonly used Teg- 
spectral type calibration by Vacca et al. (1996| ) is based 
on plane-parallel models which do not incorporate stel- 
lar winds and line blocking/blanketing. The more recent 
calibration by Martins et al. (2002| ), based on non-LTE 
line blanketed atmosphere models with stellar winds com- 
puted with the CMFGEN code of |Hillier fc Miller (1998| ), 
gives substantially lower effective temperatures than the 



Vacca et al. (1996) calibration for the same spectral type. 
We will use both calibrations. 

The result of the comparison between the observed 
spectral types and the predicted stellar properties is sum- 
marized in Table 0. The calibration by Martins et al 
(2002j)]gives effective temperatures for the observed spec- 
tral types ~ 4000 - 5000 K lower than |Vacca et al. (1996| ) . 
Even considering this most recent calibration, the CoStar 
models systematically predict a lower temperature for the 
ionizing stars than observed. A better match is found by 
using the Pauldrach SEDs, although the predicted T e g for 
G29 is still low. The opposite is true for IRAS 12073, for 
which the Pauldrach models predict a hotter star than ob- 
served. Moreover, while our photoionization models pre- 
dict that an increase in the degree of ionization is coupled 
to an increase in T e g and thus in spectral type, obser- 
vationally this trend is less clear (cf., Fig ^). In partic- 
ular, G29 and 30 Dor are both ionized by very early O 
stars, yet their ionization ratios are very different. Sim- 
ilarly, IRAS 12073 and 30 Dor are ionized by stars with 



very different spectral types but their ionization ratios are 
very similar. It is clear that the stellar spectral type as de- 
termined from optical or infrared spectroscopy is not the 
only controlling factor of the ionization structure of H n 
regions. In line with the discussion in Sect. || (cf., Fig |l]), 
we infer that the metallicity of the star and H n region is 
involved as well. At this respect, the differences in spectral 
type between G29 and 30 Dor, and in ionization structure 
between 30 Dor and IRAS 12073 may be a consequence 
of the differences in metallicity. G29 has a Ne/H abun- 
dance more than 3 times higher than 30 Dor; 30 Dor, on 
the other hand, has a lower abundance than IRAS 12073 
(cf., Fig. |). 

It is well known that a stellar wind modifies the SED 
of a star. This effect is directly coupled to the strength of 
the wind, which in turn depends on the metallicity. Metal- 
licity also controls the degree of attenuation of the UV and 
EUV stellar radiation due to line blocking and blanketing. 
Changes in metallicity can, therefore, transform the stel- 
lar spectrum to that of a cooler or hotter star than what 
is determined from the solar metallicity stellar models. 



5. Influence of the metallicity 

The ionization tracer [Ne m]/[Ne n] is compared with the 
metallicity tracer Ne/H in Fig. ||. The data show a loose 
trend of decreasing degree of ionization with metallicity. 
No selection effect is thought to influence this relation- 
ship. While some of the LMC/SMC Hn regions are among 
the brightest in these galaxies — 30 Dor (LMC) and N66 
(SMC) are powered by luminous stellar clusters contain- 
ing many early O stars — , the LMC sample also contains 
sources such as N11A, which are alike to the Hn regions 
in the Galactic sample in terms of luminosity and Lyman 
continuum flux. Besides, no selection effect is apparent in 
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the Galactic Hn regions, for which no obvious variation in ionizing stars based upon K-band spectra (Hanson et al 



the dei sity of the ionized gas or the stellar luminosity with 1996) may require spectra of stars with different metal 



Galactocentric distance is observed (cf., MHP02). There 
is some systematic effect, however, due to the influence of 
metallicity on the electron temperature of the ionized gas. 
The Galactic Ne/H abundances were derived using Bra 
to determine the H + emission associated to the nebula 
and a constant T c of 7500 K was assumed. Hence, a high 
metallicity nebula, which is expected to have a lower elec- 
tron temperature, will actually have a higher Bra emission 
coefficient than the one adopted in the evaluation of the 
Ne/H abundance, and vice versa. For reaso nable electron 



tem peratures, between 5000 and 10000 K (Shaver et al 
1983 |; 14-fflcrbach et al., 1996) ; |Dcharveng et al., 2000| ), this 



effect is less than a factor of 2 (cf., MHP02). The Mag- 
ellanic Cloud abundances, however, are not expected to 
be influenced by this systematic effect because they were 
evaluated using electron temperatures derived from opti- 
cal lines. Hence, we consider that the trend in Fig. ^ is not 
due to a selection effect, and that there is a general cor- 
relation albeit with a large spread. On average, the most 
highly ionized H n regions have the lowest metallicity, and 
vice versa. 

Giveon et al. (2001) investigated the influence of vary- 
ing the stellar and nebular metallicity on the [Ne in]/[Ne n] 
line ratio based upon a set of photoionization models com- 
puted using the Pauldrach SEDs. They found that, for a 
given T c ff, the SED hardens with decreasing metallicity 
due to the decreased line blocking in the stellar atmo- 
sphere winds. Quantitatively, for a 40 kK star, an increase 
in the nebular and stellar metallicity from 0.5Zq to 2Zq 
would reduce [Ne m]/[Ne n] by a factor of 20. We compare 
the results of these calculations to the observed correlation 
between metallicity and degree of ionization (cf., Fig. ||) 
for T eff =35, 40, 45 kK, Z=O.5Z , 1Z , 2Z Q and ?i H =800 
cm~ 3 . In view of these results, we see that the large spread 
in degree of ionization at a given metallicity observed in 
our data reflects a true variation in T e g. Moreover, the 
observed increase in the degree of ionization of the H n re- 
gions (cf., Fig. 0) reflects clearly a hardening of the SEDs 
due to the decreased metallicity and it is not necessarily 
due to a rise of the effective temperature of the ionizing 
star. 

Returning to the individual sources considered in 
Sect. |] and in view of Fig. || the Galactic Hn region 
G29, with a metallicity of roughly 2Zq, would be pow- 
ered by a star with a temperature slightly higher than 
the one indicated in Fig. ||, which are solar metallicity 
models. The Magellanic sources N66 and 30 Dor, with sub- 
solar metallicities, would be, on the other hand, ionized by 
cooler stars because of the decreased line blocking. How- 
ever, a proper comparison of the effective temperature de- 
rived from the modelling of the H n region line data with 
the effective temperature derived directly from the stellar 
spectral will require a proper calibration of stellar spectra 
for different metallicities. In particular, the 'calibration' of 



licities. Also, the comparison of the strength of hydrogen 
and helium lines with predictions of stellar atmospheres 
will require a proper treatment of the influence of metal- 
licity. 

6. The case of starburst galaxies 

Information on starburst galaxies has been collected in or- 
der to make a comparison with the H n region results. The 
data obtained for a sample of 9 galaxies observed by ISO 
are presented in Table || and are plotted in Fig. || as open 
stars. These galaxies show a similar distribution as the H II 
regions: high metallicity-low ionization, low metallicity- 
high ionization. In particular, we note that the low metal- 
licity dwarf irregular galaxy II Zw 40 coincides well with 
the Hll regions in the LMC, further reinforcing our in- 
terpretation that metallicity and degree of ionization are 
strongly interconnected. Among the low ionized galaxies 
are M 82 and NGC 253, which are the prototypical star- 
burst galaxies. NGC 3256, with a [Nem]/[Nen] similar to 
that of M 82, is a recent galaxy merger. The highly ionized 
galaxies include NGC 5253, a young, Wolf Rayet galaxy, 
and NGC 4038 /4039 (the Antennae) , where an extensive 
burst of star formation is ongoing in the overlap region 
between the two merger systems. 

The [Ne in] /[Neil] line ratio and other ionization ratios 
observed in starburst galaxies have been used to infer the 
spectral type of the ionizing stars and hence, the age of 



the starburst (e.g. Schaerer & Stasinska, 199£; Crowthci 
|ct al., 1999| ; [Thornley et al., 2000| ; |Spoon et al., 2000fl . 



However, as this discussion has emphasized, the stellar 
type is only one of the parameters influencing the ioniza- 
tion structure. Furthermore, the uncertain spectral clas- 
sification of the ionizing sources through ionization line 
ratios makes a proper age estimate of the starburst diffi- 
cult. The other factor influencing the ionization structure, 
metallicity, measures the cumulative effect of star forma- 
tion over the galaxy's history. 

7. Summary and conclusions 

In this paper, we have examined the interplay between 
the ionization structure of H n regions, the SEDs of their 
ionizing stars and metallicity. To this end, we combined 
the near- and mid-infrared ISO spectra of a set of Galac- 
tic and Magellanic Cloud Hll regions. This wavelength 
range gives us access to the [Siv] 10.5/[Slll] 18.7 /^m and 
[Nem] 15.5/[Neii] 12.8 ^tm line ratios, which are sensitive 
to the ionization structure of the nebula. The [Siv]/[Sin] 
and [Ne m]/[Ne n] line ratios are found to span a range 
of more than 2 orders of magnitude and to be highly cor- 
related. The low metallicity LMC and SMC Hll regions 
overlap the trend found for the Galactic objects at the 
high ionization end. 
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Table 3. [Neni]/[Nen] line ratios and Ne/H abundances 
for a sample of galaxies, including SgrA*. The sample is 
sorted out by increasing degree of ionization. 



Object 


Tat 1 / Ttvt 1 

[NeiliJ/[NeilJ 


Ne/H 1, 


Ref. 


name 




do -4 ) 




M83 


0.05 


2.1 


1,+ 


SgrA* 


0.06 


1.8 


2,+ 


NGC7552 


0.08 


1.8 


1,+ 


NGC253 


0.08 


1.8 


1,+ 


NGC4945* 


0.09 


1.3 


3 


M82 


0.18 


1.4 


4 


NGC3256 


0.20 


2.3 


5 


NGC 4038/39° 


0.92 


0.G 


G 


NGC5253 


4.0 


0.7 


7 


II ZW40 


12 


0.4 


1,+ 



(f) Calculated from the mid- infrared [Nell] 12.8 /im, 
[Nelll] 15.5 fim and Bra 4.0 )im lines assuming T e — 10 4 K 
and the low density limit. No extinction correction is applied. 
(J) The nature ot this galaxy is uncertain, the infrared diag- 
nostics show no evidence for the e xistence of the AGN inferred 
from hard X-ray observations (see Spoon et al., 200C). (o) ISO 



pointing centered on the overlap region between the two inter- 

actinc systems. 

REFEREN CES: (1) [Thornley et al. (2000|); (2) |Lutz et al 
(19960 (3) [Spoon et al. (20001); (4 ) |Forster Schreibcr et al 
(2001 " (5) |Rigopoulou et al. (1996| ); (6) [Kunzc et al. (1996] )~ 



(7) |Crowther et al. (1999| ); (+) ISO archive 



The dependence of the [Siv]/[Siii] and [Nem]/[Nell] 
line ratios on the T e g of the ionizing star is investigated 
using a grid of photoionization models at solar metallic- 
ity. The models are calculated using two different stellar 
atmosphere models: the CoStar SEDs from [Schaerer fc de 
Koter 1997) and the more recent SEDs from Pauldrach 



et al. (l i001| ). The interpretation of the observed ionization 



line ratios in terms of the stellar content of the H II region 
depends critically on the adopted SED. In particular, we 
note that the CoStar SEDs predict systematically much 
cooler stars than the Pauldrach SEDs for the same ioniza- 
tion conditions. As a result, practically the whole observed 
ionization range can be fit by a very narrow T c s range be- 
tween ~ 33000 and 42000 K when the CoStar models are 
adopted, while the entire range up to 55000 K is suitable 
when using the Pauldrach models. These differences orig- 
inate principally from the different treatment of the line 
blanketing/blocking and wind properties. 

For a few sources, spectral types of the ionizing sources 
are presently available from infrared (K-band) or opti- 
cal stellar spectroscopy. The comparison of these spec- 
tral types with those derived indirectly from the ionized 
gas through the [Siv]/[Siii] and [Ne III] /[Neil] line ratios 
confirms that metallicity has an important influence on 
the stellar spectra. As impetus for this, we have investi- 



gated the relation between the ionization structure and 
the metallicity comparing the observed [Neni]/[Neii] line 
ratios and Ne/H abundances. Although a large scatter is 
present, a loose correlation between the two parameters 
is found. This observed trend is compared to photoioniza- 
tion models calculated varying self-consistently the nebu- 
lar and stellar metallicity (Giveon et al., 2001). The com- 



parison shows that the observed increase in degree of ion- 
ization can be explained by the hardening of the SED due 
to a decrease of the metallicity. 

We note, moreover, that the T e g calibration of the dif- 
ferent spectral types is based on modelling the observed 
photospheric hydrogen, helium and metal lines. This cal- 
ibration is uncertain as it is very sensitive to the input 
physics of the stellar mo del used. The wide ly used cali- 
bration of early stars by Vacca et al. (1996| ) is based on 
plane-parallel, non-LTE models which do not incorporate 
stellar winds and line blanketing. It has been shown, for 
instance, that the inclusion of line blanketing, stellar wind 



and a spherical geometry in the models (Martins et al 



2002) lowers the predicted T c s substantially when com- 



pared to the calibrated by Vacca et al. (1996) for the same 
spectral type. A new calibration is needed if predictions of 
the stellar content of H II regions are to be compared with 
the observed spectral types of the ionizing stars. Such a 
recalibration will have to include the effects of metallicity. 

In view of this analysis, the main conclusions are the 
following: 

• The interpretation of the observed X +l+1 / X +l ratios 
is valid only when the stellar metallicity of the SED 
used in the modelling properly matches the metallicity 
of the local ISM. Not taking into account the effect 
of the stellar metallicity correctly can lead to wrong 
conclusions on the stellar content of the nebula and 
thus on the local stellar population. 

• The adequate treatment of the effect of stellar wind 
and line blocking and blanketing on the ionizing spec- 
tra of early O stars, and the availability of grids of 
stellar models at different metallicities, is needed. In 



this respect, the new models by Pauldrach et al. (2001 ) 
constitute an improvement over previous ones. 

• A T c ff-spectral type calibration of early stars is insuffi- 
cient if the effects of metallicity are not accounted for. 
A detailed new calibration is, therefore, necessary. 

• The determination of starburst ages in starburst galax- 
ies based on X +l+1 / X +l ratios is likely incorrect if 
metallicity is not taken into account. 

Concluding, it is important to further study the ioniza- 
tion structure of Hn regions, particularly using different 
tracers such as the Hei recombination lines. Furthermore, 
from a theoretical point of view, it may be important to 
investigate the origin of the actual spectral classification, 
such as the K band spectral classification via the C IV and 
N in photospherical lines and their dependence on the stel- 
lar metallicity. 
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